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TECHNICAL NOTE 3078 


TRANSIENT TEMPERATURES IN HEAT EXCHANGERS 
FOR SUPERSONIC BLOWDOWN TUNNELS 
By Joseph H. Judd 


SUMMARY 


In order to investigate transient temperatures in heat exchangers 
for supersonic blowdown tunnels, an analysis of the heat transfer between 
an insulated tube and the fluid moving through the tube was made for three 
inlet fluid conditions: one in which the fluid temperature was constant, 

another in which the temperature decreased exponentially, and a third in 
which the temperature of the fluid decreased linearly with time. For 
temperatures low enough to permit the neglect of radiation effects, the 
results of the analytical work give mathematical expressions whereby the 
tube and fluid temperatures may be computed as functions of axial tube 
distance and time. 

Temperatures were measured for air at a constant inlet temperature 
and at atmospheric pressure flowing through a bundle of tubes preheated 
to a constant axial temperature of 202 ° F. Excellent agreement between 
the computed and measured temperatures justifies the simplifying assump- 
tions and substantiates the analysis . 


INTRODUCTION 


One requirement for any intermittent or blowdown tunnel used in the 
st udy of jet-engine performance, in the study of heat transfer, or in tests 
of the structural integrity of airplane structures at supersonic speeds 
is the establishment of air flow at desired values of stagnation tempera- 
ture. The large amounts of heat required for this purpose are usually 
supplied by heat exc hang ers in the form of bundles of tubes or sheets 
which are preheated to approximately the desired stagnation temperature 
of the air. Although sufficient information is available for the aero- 
dynamic design of a supersonic blowdown tunnel such as that described in 
reference 1 , little information is available on the design of a heat 
exchanger. The present investigation is concerned with the transient 
heat transfer through a heat exchanger of the heat-accumulator type during 
operation of a supersonic blowdown tunnel. The effects of radiation are 
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neglected in this study and the application of the results is thus 
limited to heat exchangers in -which the temperatures are not high enough 
to cause appreciable heat losses by radiation. 

Heat exchangers of the heat-accumulator type provide a method of 
supplying very large amounts of heat for short periods of time while 
storing energy at a low rate between tests . A wide variety of materials 
have been considered, from artificially heated rocks at Peenemiinde to 
bundles of steel tubes at the preflight jet of the Langley Pilotless 
Aircraft Research Station at Wallops Island, Va. In the analysis in 
this paper, the heat is considered to be stored in bundles of tubes 
initially preheated to a constant temperature. Several air conditions 
are exam in ed by mathematical analysis. In the first condition, the air 
at a constant temperature enters a tube, as reported in reference 2, 
or a porous medium, as discussed in reference 5* I 2 - the second condi- 

tion, the air is considered to be stored in a tank under pressure . 

During the test the air from the storage tank undergoes an isentropic 
expansion before entering the heat-exchanger tubes . The third condi- 
tion represents a simplification of the preceding condition for an 
expanding gas in that a linear decrease in the temperature of the air 
entering the heat exchanger is considered. 

The basic laws of heat convection were used to derive an integro- 
differential equation which expresses the variation of the tube and air 
temperature difference as a function of time and axial distance along 
the tube. Upon Insertion of the appropriate boundary conditions for 
each case considered, expressions for the variation of the temperature 
difference were obtained. Integration of these expressions gave the 
variation of tube and fluid temperatures with respect to time and axial 
tube distance. 

In order to check the validity of the analytical method and the 
assumptions involved, an experimental investigation was made for air at 
a constant inlet temperature and at atmospheric pressure flowing thro ug h 
a tube. Two series of tests were made, one in which a single instrumented 
tube was placed in the center of a bundle of tubes and one in which the 
instrumented tube was packed in insulation. The initial temperatures of 
the tubes were 202° F and l66° F, respectively. The first series of tests 
represents the situation in a heat exchanger, whereas the second series 
represents the heating of a fluid passing through a tube where the sur- 
rounding medium is much less dense than the tube. 


SYMBOLS 


A 

c f 


internal surface area of tube per unit length, sq ft/ft 
specific heat of fluid at constant pressure, Btu/(lb)(°F) 
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Cp 

d ext 

d int 


g(t) 


h 

*0 

In 

K f 

% 

k 

l 

m 

mo 


n 

CL 

T 



T fi 


specific heat, of tube at constant pressure, Btu/(lb) (°F) 
external diameter of tube, ft 
internal diameter of tube, ft 


2»„t 

1 + Kp — — + . . . + 



r 

M 7 “\ 

(a . t r * . t) 

\*o] 

\k / 7 \k / 


7+1 


Kp 


& - «)“ 


(7 + l)(r + 2 ) . . . (7 + n) 


heat-transfer coefficient, Btu/(sq. ft)(sec)(°F) 

Bessel function of zero order 
Bessel function of order n 

ratio of heat transfer to heat capacity of fluid, l/sec 

ratio of heat transfer to heat capacity of tube, l/sec 

rate of flow through heat exchanger, Ib/sec 

length of tube, ft 

mass of fluid in tank, slugs 

initial mass of fluid in tank, slugs 

integer 

heat flux, Btu/sec 

temperature, °F 

temperature of fluid, °F 

temperature of fluid at tube entrance, °F 

temperature of external fluid before flow, °F 
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k 


Tf + 


T T 


% 

T Pt 

t 

V 

w f 

W-r, 


X 


temperature of fluid at time t, °F 

temperature of tube, °F 
initial temperature of tube, °F 
temperature of tube at time t, °F 
time, sec 

fluid velocity, ft/sec 

weight of fluid per unit length, lb/ft 

weight of tube per uni t, length, lb/ft 

axial distance along tube, ft (unless otherwise specified) 


a 



(®f f f jtepKf i(t - | n l n f KpKf - i) 


K p( t -f) +K f V 



7 ratio of specific heat of fluid at constant pressure to 

specific heat at constant volume 

0 temperature difference between tube and fluid, °F 

0 q difference between initial temperature of tube and initial 

temperature of external fluid, °F 
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0-£ difference between initial temperature of tube and temperature 

of entering fluid, °F 

p density of fluid in tank, slugs/cu ft 

Pq initial density of fluid in tank, slugs/cu ft 



HEAT TRANSFER BETWEEN A PREHEATED TUBE AND THE FLUID 


This investigation is concerned with the analysis of heat exchangers 
of the heat-accumulator type in which the heat source or sink consists of 
a bundle of tubes through which the working fluid flows. Except for the 
outer row of tubes where a small amount of heat may be transferred to 
the surrounding shell, the fluid provides the only means of heat transfer 
during operation. No heat is transferred across the outer surface of . , 
the individual tubes because al l axe at the same local temperature; that 
is, each tube is effectively perfectly insulated. 

In order to reduce the number of independent variables which must be 
considered in the analysis, the assumption was made that no radial tem- 
perature gradients exist in the tubes and that the temperature differ- 
ence between the tubes and the fluid occurs across a thin boundary layer. 
Inasmuch as the conductivity of the tubes is large compared with that of 
the boundary layer, this assumption appeaxs to be justified. For long 
tubes, turbulent flow may be considered to exist throughout the tribe. 

If the tubes have blunt ends with sharp internal lips, the flow becomes 
turbulent 3*5 diameters downstream of the entrance (ref. A) . A further 
assumption of the analysis is that the transfer of heat axially along 
the tube can be neglected. The validity of this assumption depends upon 
the temperature gradient along the tubes and upon the ratio of the tube- 
wall cross-sectional area to the internal surface area. Effectively, the 
assumption states that the heat transfer is polarized so that heat is 
transferred only to the fluid. 

For an initial tube temperature greater than the fluid temperature, 
Newton's law for the convection of heat may be written as follows: 


dq = hA0 dx 


( 1 ) 
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where q_ is the heat flux; h, the heat-transfer coefficient between the 
tubes and the fluid; A , the internal surface area of a tube; and 0, the 
temperature difference between the tubes and the fluid. Two additional 
expressions which represent the conservation of energy are obtained and 
may be expressed as follows: 


dTp 

dq - -p'p sT ^ 


( 2 ) 


dT f 

dq = w^c^ dx 

1 dt 


( 5 ) 


In these expressions Wp and Wf represent the weights of the tube and 
fluid per unit length, Cp and Cf represent the specific heats of the 
tube and fluid, and Tp and Tf represent the temperatures of the tube 

and fluid, respectively. A graphical representation of the temperatures 
and temperature derivatives is presented in figure 1. 

In order to obtain the tube and fluid temperatures , an expression 
must first be obtained for the dependent variable 0 in terms of the 
independent variables x and t. Inasmuch as 0 represents the heat 
flow potential and the fluid and tube temperatures can be expressed as 
functions of the temperature difference, this selection of the dependent 
variable was a natural one. 

The total differential for the temperature difference 0 is 


d£ = ^ dx + ST p _ ^f 
dt dx dt dt dt 


By combining equation (l) with equations (2) and ( 3 ), the following 
expressions for the temperature derivatives are obtained: 



hA 


w. 


P P°P 


( 5 ) 


hA 


0 


and 


dT f 

dt 


W-fC-f 


( 6 ) 
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In order to shorten the mathematical expressions used subsequently in 
this paper , the terms hA/w f c f and hA/v p Cp are- replaced by the coef- 
ficients Kf and K p , respectively. Physically these terms represent 

the ratios of heat transfer between the tube and the fluid to the heat 
capacities of the fluid and the tube, respectively. In order to line- 
arize the differential equations, Kf and Kp are assumed to be con- 
stants throughout the length of the tube. The partial derivative of the 
tube temperature may be obtained by differentiating equation (5) with 
respect to x and integrating with respect to t ime . Thus, 



(7) 


The values of the derivatives of the tube and fluid temperatures, 
equations (5), (6), and (7), are substituted into equation (4) in order 
to obtain an expression for 0 in terms of x and t: 


d0 

dt 



dt - V - ^ 9 


The fluid velocity V can be substituted for dx/dt. Expanding d0/dt 
gives the expression for the dependent variable 0 in terms of the 
independent variables x and t as 


V | + l = - K P V ^ t l«-(K p + Kf )e (8) 


A further condition necessary to linearize the differential equation is 
that the fluid velocity remain constant throughout the tube length. 
Rigidly then, the analysis is restricted to an incompressible fluid. 

The solution to this integrodifferential equation is the key to the 
analysis. ’Once 0 has been determined, the tube and fluid temperatures 
may be obtained from equations (5) and (6) by integration. 
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CALCULATION OF TUBE AND FLUID TEMPERATURES 
Fluid at a Constant Inlet Temperature Passing Through a Tube 
Init ially at a Constant Temperature 


BoimriFLTy conditions . - The heat-transfer condition considered in 
this section is one in which a fluid, either gas or liquid, at a con- 
stant inlet temperature flows through a tube preheated to an initial 
temperature which is constant along the length of the tube. A prac- 
tical. example of this condition is furnished in reference 1 in which 
is discussed a supersonic blowdown tunnel powered by an ejector which 
takes in air at room temperature . The boundary conditions used in the 
solution of equation (8) are determined from physical considerations of 
heat and fluid flow. 

The tube is at a constant axial temperature Tp Q before the fluid 

outside the tube, which is at a constant temperature Tf^, starts to 

flow. Before the test starts, that part of the fluid within the tube 
attains the tube temperature. Thus, the first boundary condition may 
be written 


0=0 (t < 0} 0 ^ x ^ l) (9) 


The second boundary condition is obtained by considering the tem- 
perature relationship at the tube inlet during the test. If no heat 
is transferred axially along the tube, as is assumed in the preceding 
section, a differential equation for the temperature difference at the 
inlet can be obtained. This equation is 


d0 

dt 


= -Kp0 


( 10 ) 


When t = 0, the temperature difference 0p at the tube inlet is the 
difference between T Pq and T-^. When this condition is put into the 

solution of equation (lO), the second boundary condition is determined as 


0 = 0 ± e (t ^ 0; x = 0 ) (ll) 
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Evaluation of 9 . - The integr ©differential equation for 6 as a 
function of x and t and the boundary conditions for a fluid with a 
constant inlet temperature determined in the previous sections are now 
used to evaluate 0. Equations (8), ( 9 ), and (ll) are restated for 
convenience : 

0=0 (t < Oj 0 g x ^ l) 

0 = 0 1 e"^P t (t ^ 0; x = 0 ) 


This series of equations is amenable to solution by the Laplace 
transformation. Values of 0 are given as follows: 


where 


0=0 
0 = Q ± (f> 




oQv* f 

M)J 

M*-?} 

•** f 


Equations (l2) can be shown to be a solution of the integrodifferential 
equation by direct substitution. The first boundary condition, equa- 
tion ( 9 ), is satisfied by the form of solution given by the Laplace 
transformation. When x is set equal to zero, the second boundary con- 
dition is satisfied. Values of the Bessel function are given in refer- 
ences 5 and 6. In order to provide an easy method of determining 0, 
values of f are presented in figure 2 as functions of Kpft - for 

X 

various values of Kp — . 

^ V 
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Tube and fluid temperatures . - After 9 has heen determined, the 
tube and fluid temperatures can be calculated from equations (5) and (6) . 
Although both tube and fluid temperatures may be found by integration 
of these equations , the tube temperatures are the easier ones to obtain. 
The fluid temperatures are most easily determined by subtracting 9 from 
the tube temperatures. 

If a fixed position on the tube is considered, the partial derivative 
of Tp may be changed to the total derivative and equation (5) can be 

written in the form 




9 dt 


Tp = Tp^ - Kp9qe 




L |M) 


dt 


This expression can be integrated by parts and the value of Tp then 
becomes 


T p - Tp^ - 9 pa 


(15) 


where 


a = 1 + 


(“f ff fVf -f)| v(* - f) 

K p( t_ v) +Kf 7 


As stated previously, the fluid temperature Tp is obtained by sub- 
tracting 9 (eq. (l2)) from Tp (eq. (13)). Thus, 

= t Pq - e i(“ + f) 
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Computed values of a are presented in figure 3 as functions of 
Kp (fc - for various values of K f 2 . Either equations (13) and (lk) 

and the tables of reference 6 or values of <f> and a obtained from 
figures 2 and 3 may be used to compute the temperatures fairly rapidly 
for any time and any station along the tube. 


Gas Expanded Isentropica ll y in Storage Tank and 
Flowing Through a Preheated Tube 

Boundary conditions .- In the general case of supersonic blowdown 
tunnels, an endeavor is made to keep the flow properties, pressure and 
temperature, constant in the test section during each test. When this 
condition is achieved, the gas used in the test is removed at a constant 
rate k from the tanks where it is stored under pressure. The relation- 
ship between the initial mass of gas in the t ank mQ nnfl the mass at 
any time during the test m may be written 


_m_ _ 1 kt 

“o “o 

If no heat transfer to the gas in the tank occurs, the gas in the tank 
undergoes an isentropic expansion. The mass ratio in this case is 
related to the gas temperature ratio as follows: 


1 

m_ a _p_ _ l Tf oV" 1 

“0*00 \ T fJ 

From these equations the stagnation temperature of the gas during the 
test is obtained as 



kt\ 7 - 1 
“o ) 


(15) 


Usually a valve is located between the storage tank and the heat 
exchanger to control the flow rate and to deliver the gas to the heat 
exchanger at a constant pressure. If the gas undergoes a throttling 
process in the valve, the stagnation temperature of the gas entering 
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the tubes of the heat exchanger is the same as that of the gas in the 
tank. Furthermore, the flow velocity through the heat exchanger is low 
(in order to reduce pressure losses through the exchanger) and the static 
temperature is approxim a tely the same as the stagnation temperature. 


From th i s consideration of the gas temperature at the tube inlet, 
the first boundary condition can be determined. Initially the tube has 
a constant axial temperature Tp and the gas in the tank has a tem- 


perature Tf^ . When the gas starts flowing through the tube, the gas 
temperature at the inlet Tf^ is defined by equation ( 15 ) • From equa- 
tions ( 5 ) and ( 15 ) , a differential equation is written for the temperature 
difference at the tube entrance: 


~ = -K p 0 + (7 - 1) 

at 


“o 


kb y ~ 2 

“0 / 


( 16 ) 


The value of 9 is obtained from equation (l6) and the conditions at 
the tube inlet are applied in order to obtain the first boundary condition 




(7 + l) (7 + 2) 


- *)“ 


(7 + l) (7 + 2) 


(7 + n) 


Before the test starts, the gas inside the tube attains the same 
temperature as that of the tube and the boundary condition is the same 
as that expressed by equation (9). 

Evaluation of 9 .- In order to evaluate 9, a solution must be found 
for equations (8), (9)> and ( 17 ) which are 
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V + a® „ .K P Y I* ± i 
8x 8t uq 8x 


at 


(Kp + Kf)0 


6=0 ft < 0; 0 ^ x S l) 


K "fc 

6 = [®i + T f ± s(o)] e" P - T fi s(t) (t > 0; x = 0) 

The easiest method of solution is that of the Laplace transformation. 
The solution for 0 is as follows: 



The first term of the expression for 0 ^for t - £ > 0 ) can he 

seen to he similar to that for fluid flow at a constant tube inlet tempera- 
ture, whereas the second term provides a factor which compensates for the 
variation in inlet temperature. The last part of the expression for 0 
converges fairly slowly and necessitates laborious computations involving 
very large numbers. Because 0 cannot he simply expressed as functions 

of Kp ^t - and K-p 2£ and because the computing machines available 
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do not carry a sufficient number of significant numbers, curves defining 
these functions were omitted. 

Tube and fluid temperatures .- The tube and fluid temperatures can 
be calculated by integration of equations ( 5 ) and ( 6 ), respectively. 
Inasmuch as an analytical integration seems impractical because the 
expressions are already complicated, it is advisable to deter min e the 
temperatures by graphical integration. The easiest temperature to 
determine was found to be the gas temperature. Thus, 


Tp = T. 


*i 


hA 

c^w^ 


0 dx 


0 


(19) 


The tube temperature is determined by adding 9 obtained from equa- 
tions (l 8 ) to the gas temperature obtained from equation (19)* 


A Gas With a Linear Temperature Decrease Flowing Through 

a Tube Initially at a Constant Temperature 

Boundary conditions . - In the preceding section, the problem of 
evaluating the temperatures of a tube and gas for the case where com- 
pressed gas from a storage tank flowed through a heat exchanger was 
analyzed. The variation of the stagnation temperature of the gas during 
the test is expressed by equation ( 15 ) . Use of these values produced 
expressions for Tp and T f which are difficult to evaluate. In order 

to provide a simpler method of estimating Tp and T f , the variation of 

gas temperature with time may be considered to be linear if a reasonable 
length of time is considered, that is, one-fifth or one-quarter of the 
running time allowed by the pressure energy of the stored gasj the 
variation can thus be expressed in the form 



( 20 ) 


Initially the tube has a constant axial temperature Tp Q and the 
stored gas a temperature Tf ^ . From equation ( 5 ) > an equation is obtained 
for 9 at the inlet of the tube as 


d9 

dt 


= - 



mo 


( 21 ) 
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At the start of the test the temperature difference is the constant Sq. 
When 8 is obtained from equation (2l) , the first boundary condition is 
determined so that 


8 



+ 




(t £ 0; x = 0 ) (22) 


The second bo undar y condition is determined in a preceding section 
and is expressed by equation ( 9 ) • 

Evaluation of 8 . - A flow condition exists such that a gas at con- 
stant mass rate with a linearly decreasing temperature flows through 
a tube preheated to a constant axial temperature. The problem is restated 
in mathematical terms (eqs. (8), ( 9 ), and (22)) for convenience as follows 


v S + H“- K P v / o t 5l“- ( k p + %) 9 


0=0 


(t <0; 0 <: x <: 1 ) 


0 


5 o 




^t ^ 0; x = 0 ) 

Solution of these equations for 8 was obtained by use of the Laplace 
transformation. Thus, 
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where 


P = 


1 








(Vf f) n 

(.jjn-l (-l) n - 2 l 

I 

M) + 

n! 

/ \ n-1 ‘ , 

(n - 1) IKp (n - 1) : 

n-2 

ll.’Kp 


wHT t 

i:(n - 2) IKp 



+ . • 



+ 


Although the terms in equations ( 23 ) are similar in form to those 
in equations ( 18 ) , for a gas at a variable inlet temperature the series 
in the expressions converges fairly rapidly. For this reason computation 
of 0 is relatively easy. Values of <f> are presented in figure 2 and 


values of p are presented in figure ^ as functions of 
various values of Kf — . 



Tube and gas temperatures . - Integrating equations ( 5 ) and (6) and 
using values of 0 from equation ( 23 ) gives values of tube and fluid 
temperatures . Graphical integration is mechanically the easier method 
and the temperatures of the gas are the easier to determine. 


TEST EQUIPMENT 


In order to verify the applicability of the computed temperatures 
for the fluid with constant in l et temperature, test equipment was set 
up to duplicate the conditions occurring in the heat exchanger. Ini- 
tially the tube is at a constant axial temperature greater than the 
inlet fluid temperature. Air at room temperature is drawn at constant 
velocity through the tube by a blower. This situation corresponds to an 
induction- type tunnel with an open circuit. A schematic diagram of the 
equipment is presented in figure 5- 
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The tube mounting used in these tests is illustrated in figure 6. 

The instrumented tube was a 5/18-inch hard copper tube 98 inches long 
with an inside diameter of 0.2435 inch. The interior wall of the tube 
was cleaned before instrumentation to remove oil and scale. In one 
series of tests a group of 19 tubes was used with the center tube 
instrumented (fig. 6(a)). For these tests the ins trumented tube can 
be considered to have perfect insulation because the surrounding tubes 
tend to have the same axial temperature distribution during air flow. 

The b undl e of tubes was placed in a hexagonal balsa jacket and bound 
with Paraplex impregnated Fiberglas. The diffuser and the entrance 
bell were also made of Fiberglas. In a second series of tests, a single 

instrumented tube was placed in a 2^- by 3-inch balsa insulating beam 

with an air space 7/16 inch in diameter separating the beam from the 
tube. Five equally spaced wooden washers supported the tube. A 
mahogany entrance bell and diffuser were used in this installation 
(fig. 6(b)). In these tests some heat was transferred from the insu- 
lating beam to the tube because the tube temperature decreased during 
the test. In both series of tests, the air inlet velocity was deter- 
mined from the static pressure measured at the minimum section in the 
entrance bell. 

In order to bring the tube to the desired temperature, the instru- 
mented tube, in the bundle of tubes or in the insulating beam, was 
suspended in a heating box until the equilibrium test temperature was 
reached along the tube. The heating box, 98- by 12- by 12-inch inside 
dimensions, was constructed of wood-framed celotex board. Nichrome 
wire, spaced along the box to provide a constant axial tube temperature, 
was used for the heating element. 

The air and tube temperatures were measured by thermocouples spaced 
along the tube. Thermocouple locations for the two installations are 
given in tables I and II. Iron-constantan thermocouples, 36-gage, were 
used to measure air temperatures. The thermocouple beads were 0.02 to 
0.05 inch in diameter and located along the center line of the tube. 

Tube temperatures were measured by copper-constantan thermocouples where 
30-gage constantan wire was peened into the center of the tube wall. 
Calibration of the thermocouple material for both types indicated an 
accuracy to wit hin +0.25 percent of the temperature range.' A motor- 
driven Lewis selector switch was used to switch thermocouple leads and 
transmit the individual potential to the recorders . The temperatures 
were recorded on two single-channel Brown self-balancing potentiometers, 
one for air temperatures and the other for tube temperatures. The basic 
accuracy of these ins tr ume nts was "tO.25 percent of full-scale deflection. 
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TESTS 


In preparation for each test, the tube was placed in the heating 
box until the test temperature was obtained. The uniformity and magni- 
tude of the temperature were checked by observing the tube thermocouple 
readings. Dur ing this period, the entrance bell was closed. Immediately 
before each test the duct was opened and the recorders started. Then 
the blower was started and the air valve set to obtain the desired inlet 
velocity. Air at roam temperature was drawn through the tube and heated. 
This inlet air remained at essentially a constant temperature. The test 
procedure was the same for the single insulated tube and the bundle of 
tubes . 

The static pressure in the entrance bell was measured on an electri- 
cally driven micromanometer. The thermocouple potential for the tube and 
fluid were recorded by self-balancing potentiometers. The barometric 
pressure was used as the value of total pressure during a test. 

Data are presented for a bundle of tubes at an inlet velocity of 
75 feet per second. The initial tube temperature was 202° F and the 
initial air temperature was 86° F. The inlet Reynolds number, based on 
the internal diameter of the tube of 0.2^35 inch, was 8,620 for these 
tests . 

For a single tube in insulation, v a test was made at an inlet velocity 
of 103 feet per second and a Reynolds number of 11,780. The initial tube 
temperature was l66° F and the entering air temperature was 74° F. 


ACCURACY 


The accuracy of the individual components of the temperature measure- 
ment system is given in the section entitled "Test Equipment." Calibra- 
tions indicated an error of to . 25 percent of the temperature range for the 
thermocouple material and t 0.25 percent of full-scale deflection for the 
recorders. Thus, mec hani c ally the temperature error totals iO .50 percent 
of f ull -scale deflection. However, a further source of error occurs in 
printing and reading the records. This error is estimated to be to. 5 0 F. 
Thus, the experimental error involved in the temperature measurements 
varies from 12.70 percent for the low temperatures to I 1 . 0 O percent for 
the high temperatures. 


RESULTS ADD DISCUSSION 


In prece ding sections, the heat transfer between a fluid and a tube, 
initially preheated to a constant axial temperature, was analyzed and the 
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temperature difference 9 was determined for a fluid with, several types 
of temperature variations at the tube inlet. In addition to being appli- 
cable to bundles of tubes, the results may, when the appropriate heat- 
transfer coefficients are chosen, be applied to heating elements which 
form a continuous surface parallel to fluid flow such as plates. Fur- 
thermore, large numbers of staggered rods normal to the flow, porous 
material, and similar heat-accumulator material may also be evaluated. 

The tube and fluid temperatures were obtained by integration of 0, 
analytically in the case of a constant fluid inlet temperature and 
graphically for linearly and isentropically decreasing gas temperatures . 

In the analysis, heat losses due to radiation were neglected and the 
results are limited to those cases in which the heat -exchanger tempera- 
tures are low enough to justify this condition. 

Correlation of computed and experimental temperatures .- In the 
analysis of the tube and, fluid temperatures, one of the major assumptions 
is that the tube and fluid properties are invariant with time and dis- 
tance along the tube, that is, independent of temperature variation. It 
was shown in reference 4 that turbulent flow starts 3-5 diameters from 
the inlet of blunt-edged tubes. Thus, the use of an average turbulent 
heat-transfer coefficient (ref. 7 ) over the entire length of the heat- 
exchanger tubes is justified in evaluating the coefficients Kp and Kf . 

From the definitions of these terms, it can be seen that, for all prac- 
tical purposes, Kp is constant along the tube. However, Kf is a 
function of the fluid density; for liquids Kf can be considered constant, 
but for gases a large variation occurs. Two methods of using the fluid 
properties were tried in computing the coefficients. In one, the prop- 
erties of air at the tube inlet temperature were used, and in the other, 
the properties of air were determined for an average temperature between 
the tube temperature and the air inlet temperature. Of these methods of 
determining coefficients, use of the average temperature provided closest 
agreement between experimental and computed temperatures . Inasmuch as 
the correlation was good over the entire length of the tube and over a 
fairly long time interval, the assumption that Kf was constant is 
justified. 

In order to check the validity of the analytical method and the 
assumptions involved, tests were made in which air at a constant entrance 
temperature was drawn through a tube initially preheated to a constant 
axial temperature. This condition was tested because of the simplicity 
of the test equipmient required and because, in dealing with a gas, the 
greatest departure from analytical conditions prevailed. 

Multiple-tube bundle ,- A comparison of computed and measured tube 
temperatures is presented in figure 7(a) for an instrumented tube in 
the center of a bundle of tubes . The computed temperatures presented 
are based on the properties of air evaluated at the average air 
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temperature . Essentially, the measured temperatures are those of a 
perfectly insulated tube as required by the analysis . . The computed tube 
temperatures are very close to the experimental values over the entire 
length of the -tube and for the time interval shown. Such agreement 
seems to justify f ully the analysis presented in the preceding sections 
of this paper. 

The values of air temperatures given by equation (l 4 ) represent the 
mean air temperature across the tube. The air thermocouples were located 
along the center line of the tube and, thus, give the minimum temperature 
at their respective axial stations. For the purpose of correlating the 
computed air temperatures, a l/7-power temperature distribution was taken 
across the tube because the flow was considered to be turbulent. A rela- 
tionship between the mean temperature and the minimum air temperature was 
obtained by integration of the air temperature distribution across the 
tube. The values of minimum air temperatures at the center of the tube 
and the mean air temperature are plotted in figure 70 3 ) • The computed 
and measured temperatures at the center of the tube can be seen to agree 
except near the tube inlet where a fully developed turbulent temperature 
distribution does not occur. However, it is beyond the scope of this 
investigation to determine the temperature distribution across the tube 
at these stations. 

Single insulated tube .- In order to show the magnitude of error 
involved when heat transfer resulting from imperfect insulation exists , 
data are presented in figure 8 for a single tube in insulation. The 
computed tube and air temperatures presented in figures 8(a) and 8(b) 
are based on the properties of air at the average temperature. Both 
tube and air temperatures had higher experimental values than corre- 
sponding computed values although fair correlation was obtained at the 
downstream end of the tube. Inasmuch as the heat content of the insula- 
tion was not attributed to the tube because of the nonuniformity of the 
heat flow path, this condition can represent an appreciable error. If 
a uniform insulation about the tube were used, inclusion of the mass 
of the insulating material in the coefficient Kp should lower the dif- 
ferences between the computed and the experimental values . 


CONCLUDING REMARKS 


An analytical method has been presented for the computation of tube 
and fluid temperatures in a heat exchanger consisting of bundles of tubes 
preheated to a constant axial temperature. Results have been obtained 
for three cates of interest in the design of supersonic blowdown tunnels: 
the cases in which the air inlet temperature is constant, decreases 
exponentially, and decreases linearly with time of test. Although the 
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conditions of the paper deal with heating, the results can be applied 
equally well to cooling processes. 

Computed and experimental temperatures for air at a constant 
entrance temperature and at atmospheric pressure flowing through a 
bundle of tubes preheated to 202° F have shown excellent correlation. 
This agreement justifies the method of temperature calculation. 

For the case of air at a constant entrance tenperature anrl at 
atmospheric pressure flowing through a single insulated tube preheated 
to 166° F, the experimental values of temperatures were hi ghp-r than 
computed values. However, at the tube exit, the differences were not 
too great for practical use. 


Langley Aeronautical Laboratory, 

national Advisory Committee for Aeronautics, 
Langley Field, Va., January 28, 195^. 
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TABLE I 

LOCATION OF THERMOCOUPLE IN TUBE BUNDLE 


Location of thermocouple 
for measuring temperature 
of air, x, in. 

Location of thermocouple 
for measuring temperature 
of tube, x, in. 

0 

0.125 

2 

2.125 

4 

4.125 

8 

8.125 

16 

16.125 

26 

26.125 

36 

36.125 

56 

56.125 

76 

76.125 

96 

95.875 
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TABLE II 

LOCATION OF THERMOCOUPLE IN SINGLE INSULATED TUBE 


Location of thermocouple 
for measuring temperature 
of air, x, in. 

Location of thermocouple 
for measuring temperature 
of tube, x, in. 

0 

0.125 

4 

4.125 

10 

10.125 

18 

18.125 

28 

28.125 

38 

58.125 

48 

48.125 

58 

58.125 

68 

68.125 

78 

78.125 

96 

95.875 
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Figure 5 .- Diagram of experimental setup for single tube with air flow 

at constant inlet temperature. $ 
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(a) Multiple- tube installation for air flow at constant inlet temperature . 


Figure 6.- Illustration of tube mounting. 
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(ti) Single-tube installation for air flow at constant inlet temperature. 

Figure 6.- Concluded. 
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(a) Tube temperatures . 

Figure 7*- Comparison of computed and measured temperatures for multiple' 
•tube bundle. (All values of x in inches.) 
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(t>) Air temperatures. 
Figure 8.- Concluded. 
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